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ABSTRACT 

Direct imaging of exoplanets requires both high contrast and high spatial resolution. Here, we 
present the first scientific results obtained with the newly commissioned Apodizing Phase Plate coro- 
nagraph (APP) on VLT/NACO. We detected the exoplanet /3 Pictoris b in the narrow band filter 
centered at 4.05 /im (NB4.05). The position angle (209.13°±2.I2°) and the projected separation to its 
host star (0."354 ± 0".0I2, i.e., 6.8 ± 0.2 AU at a distance of 19.3 pc) are in good agreement with the 
recently presented data from Lagrange et al. (2010). Comparing the observed NB4.05 magnitude of 
11.20 ± 0.23 mag to theoretical atmospheric models we find a best fit with a 7-10 Mj up i ter object for 
an age of 12 Myr, again in agreement with previous estimates. Combining our results with published 
V photometry we can compare the planet's [V — NB4.05] color to that of cool field dwarfs of higher 
surface gravity suggesting an effective temperature of ^1700 K. The best fit theoretical model pre- 
dicts an effective temperature of ^1470 K, but this difference is not significant given our photometric 
uncertainties. Our results demonstrate the potential of NACO/ APP for future planet searches and 
provides independent confirmation as well as complementary data for j3 Pic b. 

Subject headings: stars: pre-main sequence — stars: formation — planets and satellites: formation — 
planets and satellites: detection 
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1 Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, un- 
der program number 060.A-9800(J). 



1. INTRODUCTION 

The first images of exoplanets around stars were pub- 
lished in the last two years (HR8799 b, c, d - Marois et al. 
(2008); Fomalhaut b - Kalas et al. (2008); 1RXS J1609- 
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2105 - Lafreniere et al. (2008, 2010)) following the dis- 
covery of a planetary mass companion to the brown dwarf 
2MASS1207-3932 in 2005 (Chauvin et al. 2005). For the 
planets around HR8799 follow-up observations at 3-5/xm 
revealed evidence for non-equilibrium chemistry in the 
planetary atmospheres (Janson et al. 2010; Hinz et al. 
2010) highlighting the potential of multi-wavelength im- 
agery of extrasolar planetary systems to constrain phys- 
ical properties in comparison to theoretical atmospheric 
models. 

Lagrange et al. (2009a) detected a planetary mass can- 
didate around the young A-type star ft Pictoris (A5V, 
12±| Myr, 19.3±0.2 pc, 1.75 M Q ; Zuckerman et al. 2001; 
Crifo et al. 1997) in L' images from 2003. After a non- 
detection based on data obtained in early 2009 (Lagrange 
et al. 2009b), Lagrange et al. (2010) recently re-detected 
the planet ft Pictoris b which had moved on the other 
side of the star. The authors constrained the semi-major 
axis to 8-15 AU (the smallest of all imaged exoplanets 
today) and derived a mass estimate of 9±3 Mj up it r- 

Here, using the newly commissioned apodizing phase 
plate (APP) on VLT/NACO (Kenworthy et al. 2010; 
Codona et al. 2006), we confirm the detection of ft Pic- 
toris b with an independent data set and provide com- 
plementary photometry. 

2. OBSERVATIONS AND DATA REDUCTION 

The data were obtained on 2010-04-03, during the com- 
missioning of the APP with the AO high-resolution cam- 
era NACO (Lenzen et al. 2003; Rousset et al. 2003) 
mounted on ESO's UT4 on Paranal. The APP is de- 
signed to work in the 3-5 ^tm range where it enhances 
the contrast capabilities between ~2-7 X/D on one side 
of the PSF. We used the L27 camera (~ 27 mas pixel -1 ) 
with the visible wavefront sensor. All images were taken 
in pupil stabilized mode with the NB4.05 filter (A c = 
4.05 ^m, AA = 0.02 ^m) 2 . We used the "cube mode" 
where all data frames are saved individually. To ensure 
that no frames were lost we only read out a 512x512 
pixel sub-array of the detector. As we knew the position 
angle where the planet was to be expected (Lagrange 
et al. 2009a, b) we rotated the camera by 140° so that 
the planet would appear in the "clean", high-contrast 
side of the APP PSF (Fig. 1, left panel). 

We obtained six data cubes of ft Pictoris and di- 
rectly thereafter six data cubes of the PSF reference 
star HR2435. HR2435 has already been used in previous 
studies as reference star for ft Pictoris (e.g., Lagrange 
et al. 2009a) as one can obtain data at matching paral- 
lactic angles. Also, the star is close in the sky and has a 
similar brightness providing comparable AO correction. 
For both sources, each data cube was taken at a slightly 
different dither position following a 3-point dither pat- 
tern which was repeated twice. The on-source integra- 
tion time was 20 minutes each. Table 1 summarizes the 
observations and the observing conditions. We chose to 
saturate the core of the stellar PSFs but we note that the 
APP reduces the peak flux in the PSF core by roughly 
40% so that the comparatively long exposure time of 1 
sec did not lead to saturation effects outside of the inner 

2 The APP has been optimized to work with the NB4.05 and 
IB4.05 filter but it can also be used with the L' filter (Kenworthy 
et al. 2010). 



~5 pixels. 

For the photometric calibration we also obtained un- 
saturated images of ft Pictoris (~50% Full Well) prior 
to the science observations described above. We used 
the same observing strategy but decreased the detector 
integration time to 0.2 sec. 

The data reduction was done using self-developed IDL 
routines. The following steps were applied to all three 
datasets (i.e., ft Pictoris unsaturated and saturated im- 
ages, HR2435 saturated images). First, in order to elim- 
inate the sky background emission, we subtracted from 
each individual frame the corresponding frame from an- 
other cube taken at a different dither position (e.g., Cube 
2 frame 10 - Cube 3 frame 10) and vice versa. As the 
first two frames in each cube always showed a higher 
detector noise level they were disregarded. We also dis- 
regarded frames where the AO correction was poor (flux 
measured in the PSF core less than 50% compared to 
the previous frame). Then, bad pixels and cosmic ray 
hits deviating by more than 3-cr in a 5x5 pixels box 
were replaced with the mean of the surrounding pixels. 
After this we continued as follows: The unsaturated ft 
Pictoris images were aligned and cube- wise median com- 
bined, yielding six final images on which we performed 
photometry (Fig. 1, left panel). The alignment of the 
images was done using cross-correlation where the opti- 
mal shift between two images was determined with an 
accuracy of 0.1 pixel. The same procedure was applied 
to the saturated HR2435 images yielding six individual 
reference PSF images (i.e., one per cube; Fig. 1, mid- 
dle panel). The individual frames were not de- rotated 
to the same field orientation before the combination, so 
that all static telescope aberrations remained constant 
throughout the cubes and in the final images. For the 
saturated frames of ft Pictoris we determined the paral- 
lactic angle for each individual frame by linear interpola- 
tion between the first frame and the final frame in each 
cube 3 . Using cross-correlation, all frames from all cubes 
were then aligned to the same reference image for which 
we used the final, median-combined image of the first 
cube of the unsaturated ft Pictoris exposures. There- 
after we aligned, scaled and subtracted one of the six 
final HR2435 reference PSFs from each individual satu- 
rated ft Pictoris frame. The choice of the reference PSF 
was a trade-off between matching parallactic angle and 
observing conditions in the individual cubes. The best 
results (i.e., strongest signal of the planet, least residu- 
als) were obtained using reference PSF 2 for the frames 
in cubes one and two, PSF 3 for cubes three and four, 
and PSF 5 for cubes five and six. The scaling was done in 
the bright, righthand side of the PSF where the diffrac- 
tion rings were clearly visible. We scaled the reference 
PSF by minimizing the mean in the subtracted images in 
a semi-annulus centered on the star with an inner radius 
of 15 pixels, an outer radius of 27 pixels (i.e., covering the 
3 rd and 4 th Airy ring) and an opening angle of 150°. Fi- 
nally, we dcrotatcd all PSF-subtractcd frames to match 
the parallactic angle of the first frame and created a me- 
dian combined final image. 

3 Per default, only the parallactic angle at the beginning and at 
the end of each cube are saved in the fits header in NACO's cube 
mode. 
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3. RESULTS 

3.1. The detection of (3 Pic b 

The right-hand panel in Fig. 1 shows the final image 
of the data reduction process. The exoplanet is clearly 
detected in the left side of the image. We conducted sev- 
eral tests to exclude the possibility of a false detection: 
(1) We created final images for each individual cube by 
derotating and combining the respective PSF subtracted 
images. The planet was detected in the final images of all 
six cubes. (2) Comparing the position angle of the exo- 
planet in the final images of the individual cubes revealed 
the expected rotation introduced by the pupil tracking 
mode. (3) Speckles can appear along the spider arms 
holding the secondary mirror of the telescope, but the 
nearest arm of the telescope spider was > 20° away from 
the position of the planet. (4) We implemented a second 
independent data reduction pipeline based on the LOCI 
algorithm (Lafreniere et al. 2007). As the small field ro- 
tation in our data set did not allow us to construct a 
reference PSF directly from the (3 Pic images (Angular 
Differential Imaging; e.g., Marois et al. 2006) we used the 
HR2435 data set. For each (3 Pic frame we constructed a 
reference PSF based on a linear combination of the final 
HR2435 images so that the residuals in the high-contrast 
side of the PSF were minimized. With this approach we 
also recovered the signal of the exoplanet at the same 
position. 

The projected separation between star and planet is 
0."354± 0".012 (i.e., 6.8 ± 0.2 AU at a distance of 19.3 
pc) and the position angle is 209.13° ± 2.12° (East of 
North). The values are mean values derived from the six 
final PSF subtracted cube images (see, test (2) above) 
and the errors include the corresponding standard devi- 
ation of the mean and systematic uncertainties from the 
alignment of the images. The position of the star in the 
unsaturated reference image used for the alignment and 
the position of the planet in the six final images were 
determined by fitting a 2-dimensional Gaussian to the 
respective source. For our final astrometric numbers we 
had to rely on the plate scale and field orientation pro- 
vided in the header of the images as we did not observe 
an astrometric calibrator. However, potential deviations 
are expected to be small compared to our errors (La- 
grange et al. 2009a). Comparing the final numbers to 
the results from Lagrange et al. (2010) we find that the 
position angle is in very good agreement and the sepa- 
ration appears to have increased as expected from the 
planet's orbital motion. However, the short time base- 
line between our data and the data from Lagrange et al. 
(2010) does not allow us to put any new constraints on 
the planet's actual orbit. 

3.2. Photometry and color of (3 Pic b 

To derive the relative photometry between the planet 
and its star we used an aperture with a radius of 2 pixels 
and computed the mean flux of (3 Pictoris in the six final 
unsaturated images and the mean flux of the planet in 
the six final PSF subtracted cube images. The flux de- 
rived from the unsaturated images was scaled to account 
for the difference in the integration time. Each flux mea- 
surement was corrected for residual background emission 
by measuring and subtracting the mean flux per pixel in 
a semi-annulus centered on the planet as well as on the 



star. The semi-annulus covered only the high-contrast 
side of the objects and had an inner and outer radius of 
5 and 8 pixel for the planet and of 18 and 25 pixel for 
the star. The contrast between the star and the planet in 
the NB4.05 filter amounts to Am 4 .o5 = 7.75 ± 0.23 mag. 
The error was derived from error propagation taking into 
account the standard deviations of the mean of both the 
flux of the planet and the flux of the unsaturated (3 Pic- 
toris images. 

In the L' filter the contrast between the planet and 
the star is Ami/ = 7.7 ± 0.3 mag 4 which, combined 
with the stars apparent magnitude of L' = 3.454 ± 0.003 
mag (Bouchet et al. 1991), translates into an appar- 
ent magnitude for the exoplanet of my = 11.15 ±0.3 
mag. Since an A5V star has an \L' — M) color of ^0.01 
mag we can assume that the intrinsic \L' — NB4.05] 
color of j3 Pic is negligible. Based on the observed con- 
trast in the NB4.05 filter the exoplanet's apparent mag- 
nitude is then mpjB4 05 = 11.20 ± 0.23 and its color 
[V - NB4.05] = -0.05 ± 0.38 mag. 

3.3. Color comparison to cool field dwarfs 

Cushing et al. (2008) published NIR spectra for cool 
L and T field dwarfs that cover the NB4.05 filter but 
terminate before the long wavelengths cut-off of ESO's 
L' filter. However, for some of these objects Leggett 
et al. (2002) published L' photometry. To obtain the 
magnitudes in the NB4.05 filter we used the filter trans- 
mission curve and convolved it with the spectra from 
Cushing et al. (2008). The zero point of the filter was 
derived from a Kurucz model of an A0V star. Using the 
L' filter transmission curve we first scaled the Kurucz 
model to the observed flux density of Vega in the L' fil- 
ter. For this we assumed m^f sa ' = mj* = 0.03 mag 
and Fj c s a = 246.105 Jy or F^ cga = 5.219 • 10~ n W wT 2 
fim^ 1 in the ESO V filter. Then we measured the flux 
density of the A0V star in the NB4.05 filter and derived 
the zero point assuming the star has the same magnitude 
here. The derived flux densities of the field objects could 
then be converted into magnitudes and we determined 
the [L 1 — NB4.05] color for the objects. The error in the 
magnitude in the NB4.05 filter was derived from com- 
puting a minimum and a maximum magnitude by either 
subtracting or adding the mean error of the spectra in 
that wavelength range. Whatever resulted in a bigger 
deviation from the mean magnitude was used as (con- 
servative) error bar. Finally, as the V photometry from 
Leggett et al. (2002) was done in the MKO photometric 
system and not in the ESO system we applied a first or- 
der correction by subtracting the zero point offset of 0.02 
mag 5 from the quoted photometric measurements. 

In Fig. 2 we plot the spectral type of the field objects 
against their color (black points). Data were available 
for spectral types LI, L3.5, L7.5, T2 and T5 6 . Fitting 
a straight line to the data points we find that the color 
of f3 Pic b is most similar to that of a field object with 

4 This figure was obtained from the October 2009 dataset in 
Lagrange et al. (2010) as well as from the November 2003 dataset 
in Lagrange et al. (2009a). 

5 see, http://www.gcmini.edu/sciops/instrumcnts/midir- 
rcsourccs/imaging-calibrations/fluxmagnitude- conversion 

6 Optical spectral types differ slightly from NIR spectral types 
for 3 sources. However, the outcome of this analysis remains un- 
changed if we use the NIR spectral types instead. 
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a spectral type of L4. However, spectral types from LO 
to L7 are consistent with the observed color given the 
large error bars. For field objects around spectral type 
L4, Cushing et al. (2008) found effective temperatures 
of 1700 K which is the same temperature one obtains 
using the spectral type - temperature relation derived 
by Golimowski et al. (2004). We note, however, that the 
surface gravity of the young exoplanet is expected to be 
much lower than that of the field objects (log g « 4.0 
rather than 4.5-5.5 for field L/T dwarfs; Cushing et al. 
2008). 

3.4. Comparison to atmospheric/evolutionary models 

In Fig. 3 we compare the observed magnitudes to the- 
oretical isochrones for low-mass objects based on the 
DUSTY models from Chabrier et al. (2000). These evo- 
lutionary models are based on a spherical collapse model 
where the objects start initially with very large radii and 
their internal energy is dominated by contraction. Us- 
ing a distance of 19.3 pc and assuming an age of 12 Myr 
we find a mass of 7-10 Mj up it e r which agrees with the 
6-12 Mj U pitcr derived from the L' filter. Here, the mass 
range is only determined by the photometric uncertain- 
ties. The most likely mass based on the NB4.05 magni- 
tude is 8 Mj U piter which corresponds to T e s w 1470 K 
and log g w 4.0. Assuming an age of 20 Myr this mass 
would increase to 11 Mj up it r with T e s rss 1425 K and log 
g w 4.1, while for 8 Myr we would estimate 6 Mj up i tcr 
with T off w 1380 K and log g w 3.9. 

3.5. Contrast curve analysis 

Using the final image (right panel, Fig. 1) we computed 
a contrast curve showing the detection limit for potential 
additional companions as a function of radius (Fig. 4). 
To derive the curve we computed the standard deviation 
a of the pixel values in semi-annuli (5 pixel width) as a 
function of radial distance from the center in the high- 
contrast side of the image. The signal and the noise of a 
hypothetical companion can be written as 

S = n-r]-<7 (1) 

and 

NRiy/n~-o- (2) 

with n being the number of pixels in the chosen photo- 
metric aperture (here: n — ir ■ 2 2 w 12.56) and rj being 
a factor that depends of the chosen S/N ratio one wants 
to achieve. For a S/N of 5, 7/ can be derived from 

S 

— =5~Vn-V (3) 

The ratio between the signal of a planet with a S/N of 
5 (equation (1)) and the mean flux of beta Pictoris com- 
puted in a 2 pixel aperture relates to a difference in mag- 
nitude which we compute as a function of cr(r) (Fig. 4). 



As the noise of the residuals is not perfectly Gaussian 
in all annuli (see also, e.g., Kasper et al. 2007) we chose 
to plot the conservative 5-cr limit only. A more sophis- 
ticated analysis, i.e., the insertion of fake planets, may 
result in a somewhat different contrast curve. 

4. DISCUSSION 

j3 Pictoris b is an exoplanet whose estimated mass, age 
and separation are in agreement with predictions from 
core accretion planetary formation models (see, Lagrange 
et al. 2010, and references therein). It is interesting to 
note that dynamical studies analyzing the sub-structure 
of the remnant debris disk around j3 Pic predicted the 
mass and the orbit of the planetary companion before it 
was imaged (e.g., Mouillet et al. 1997; Augereau et al. 
2001). It was shown that an object with a mass sig- 
nificantly higher than ~10 Mj up i tcr would create much 
larger disk asymmetries. Thus, the dynamical analysis 
provides further support for the mass estimates derived 
from the observed photometry and the evolutionary mod- 
els. This in turn led Lagrange et al. (2010) to the conclu- 
sion that the "cold start" evolutionary models of young, 
giant planets (Fortney et al. 2008) are not in agreement 
with (3 Pic b. These models predict the planet to be 
significantly less luminous for its age and dynamically 
predicted mass. Our new flux point at 4.05 ^m confirms 
the brightness of the source and supports the argument 
from Lagrange et al. (2010). However, given the com- 
paratively large uncertainties in the observed fluxes we 
refrain from undertaking a more comprehensive compar- 
ison to additional sets of evolutionary and atmospheric 
models (e.g., Burrows et al. 2006). For instance for an 
assumed age of 12 Myr the models from Chabrier et al. 
(2000) predict [I/-NB4.05] < 0.3 mag for all objects be- 
tween 4 and 100 Mj up it cr . Also, the apparent difference 
in effective temperature of f3 Pic b derived in sections 3.3. 
and 3.4. is not statistically significant. Given the error 
bars in Fig. 2, T e s as high as ^1700 K is still consistent 
with the observed magnitudes within 2-a. 

Future observations at different wavelengths will help 
us to better constrain the atmospheric parameters and 
composition of f3 Pic b and to check whether the prelimi- 
nary result, that the models predict a lower T e g than the 
comparison to field dwarfs, persists. A frequent monitor- 
ing of the exoplanet 's position will eventually allow us to 
determine the orbital parameters with higher accuracy. 

This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France. We are indebted 
to Michael Cushing, Isabelle Baraffe, Udo Wehmeier and 
the ESO staff on Paranal, in particular Jared O'Neal, 
for their support. We also thank the referee for useful 
comments and suggestions. 

Facilities: VLT:Yepun (NACO) 
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TABLE 1 

Summary of observations on 2010-04-03 

Parameter fS Pictoris HR2435 

UT start 00:39:31.1 01:34:21.3 

NDIT x DIT a 200 x 1 s 200 x 1 s 

NINT* 6 6 

Parallactic angle start 69.469° 69.451° 

Parallactic angle end 74.492° 74.499° 

Airmass 1.38-1.44 1.42-1.48 

Typical DIMM seeing 0.70-0.90" 0.60-0.75" 

(EC) C 29.5-44.4 % 22.7-42.9 % 

(r ) d 4.3-7.9 ms 4.0-7.8 ms 



a NDIT = Number of detector integration times (i.e., 
number of individual frames); DIT = Detector inte- 
gration time (i.e., single frame exposure time) 

b NINT = Number of integrations or data cubes (2 
X 3 dither positions = 6 data cubes) 

c Average value of the coherent energy of the PSF per 
data cube. Calculated by the Real Time Computer 
of the AO system. 

d Average value of the coherence time of the atmo- 
sphere per data cube. Calculated by the Real Time 
Computer of the AO system. 



NACO/APP 4-micron observations of j3 Pictoris b 

Logarithmic scale Logarithmic scale 



Linear scale 
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Fig. 1. — Left: Median combination of one cube of unsaturated exposures of Pictoris used to determine the photometry. The effect of 
the APP can be seen in the left side of the PSF where the diffraction rings arc effectively suppressed increasing the contrast between 2 and 
7 X/D. Middle: Median combination of one cube of saturated exposures of the PSF reference star HR2435. Right: Median combination 
of all PSF subtracted science exposures, f) Pictoris b is indicated by the arrow. The right side suffers from subtraction residuals as does 
the central region of the PSF which has been masked out. 
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Color vs. spectral type 
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FlG. 2. — The [L' - NB4.05] color as function of spectral type derived from cool field dwarf (black points). The dashed line is a linear fit 
to the field dwarf data. The position of f) Pic b on the fitted correlation is shown in red. The error bars are the root-sum-squares of the 
photometric errors in the individual filters. The shaded area covers the range of spectral types consistent with the color of the exoplanet 
given the error bars. 
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Comparison with evolutionary models 
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Fig. 3. — The apparent L' and NB4.05 magnitudes of /3 Pic b (black cross) compared to theoretical isochrones based on the DUSTY 
models from Chabrier et al. (2000). Circles, triangles and boxes denote objects with 5, 10, and 15 Mj up it C r, respectively. The age of the 
objects is color-coded with red symbols being 8 Myr, black being 12 Myr, and blue being 20 Myr. The dashed black line shows the 12 Myr 
isochrone. 
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Contrast curve 
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Fig. 4. — 5-<r Contrast curve derived from the high-contrast side in final PSF-subtracted image shown in the right panel of Fig. 1. The 
dashed horizontal lines correspond to masses as predicted by the DUSTY models (Chabrier et al. 2000) for an assumed age of 12 Myr. The 
red cross shows the position of /3 Pic b. 



